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Abstract 

We exhibit an interaction-driven metal-insulator quantum phase transition in a holo- 
graphic model. Use of a helical lattice enables us to break translation invariance while 
preserving homogeneity. The metallic phase is characterized by a sharp Drude peak 
and a d.c. resistivity that increases with temperature. In the insulating phase the Drude 
spectral weight is transferred into a 'mid-infrared' peak and to energy scales of order the 
chemical potential. The d.c. resistivity now decreases with temperature. In the metallic 
phase, operators breaking translation invariance are irrelevant at low energy scales. In 
the insulating phase, translation symmetry breaking effects are present at low energies. 
We find the near horizon extremal geometry that captures the insulating physics. 



Metal-insulator transitions driven by electron interactions present serious theoretical 
challenges and yet are a key feature in the phase diagram of several families of materials, 

Metal- 

including the cuprate superconductors [1, 2]. As the system transitions to an insulating 

INSULATOR 

state, the d.c. conductivity drops dramatically, while in the optical conductivity the Drude 

TRANSITIONS 

peak gives way to a gap in the low energy spectral weight [1, 3]. In the vicinity of the metal- 
insulator transition one can encounter 'bad metals'. These are characterized by metallic 
behavior in the absence of a coherent Drude peak and have resistivities exceeding the Mott- 
Ioffe-Regel saturation limit [4, 5, 6]. Figure 1 shows a cartoon of the optical conductivity in 
conventional metals, bad metals and Mott insulators. It has been argued that bad metals 
cannot admit a quasiparticle description because such quasiparticles would necessarily have 
a mean free path shorter than their Compton wavelength [7]. In some cases at least, 
the reorganization of the charge-carrying degrees of freedom, from itinerant to localized 
excitations, across an interaction-driven metal-insulator transition is likely to involve strong 
coupling dynamics that is inaccessible to conventional theoretical tools. 
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Figure 1: Progressive loss of Drude coherence in the optical conductivity. From top to 
bottom on the left: conventional metals, bad metals and Mott insulators. 



The holographic correspondence allows a description of certain strongly interacting me- 
dia that eschews quasiparticles from the outset [8] . This approach has been demonstrated to 
offer new qualitative and quantitative insights into transport in quantum critical particle- Holography 
hole symmetric theories, e.g. [9]. In this paper we will shall use holographic duality to 
give an inherently strongly interacting description of a metal-insulator transition in a finite 
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density theory. Our results for the optical (in figure 4) and d.c. (in figure 5) conductivities 
capture key qualitative features of the experimental data on metal-insulator transitions, 
including bad metallic regimes [4, 5, 6] and spectral weight transfer [3, 4, 10, 11]. 

The primary obstacle to the realization of a genuine insulating phase in holography has 
been the fact that such a phase must break translation invariance. In metallic phases with a 
sharp Drude peak, the effects of momentum non-conservation are an irrelevant deformation 
of the continuum theory at low energies that may be treated perturbatively [12]. Thus the 
effects of irrelevant disorder [13, 14] and lattices [12] on charge transport at low energies and 
temperatures are easily incorporated into holographic models. The modeling of bad metals 
and insulators, however, requires the breaking of translation invariance to be treated non- 
perturbatively at the lowest energy scales. Recent impressive steps in this direction have 
been taken in [15, 16], where a lattice was fully incorporated into holographic models. 
However, the regimes studied thus far have all exhibited a conventional Drude peak. 

An important technical innovation of our work will be to break translation invariance 
while retaining homogeneity of the system. 1 This will enable our holographic study of local- 

Helical 

ization physics to operate at the level of ODEs in the bulk, while remaining nonperturbative 

LATTICES 

in the lattice strength, as opposed to the technically involved PDEs that are typically nec- 
essary once boundary translation invariance is broken. Specifically, we turn on the following 
boundary background fields (couplings) 

AW=ndt, B (0) = Aw 2 . (1) 

Here fi is the chemical potential for the electric charge while is a source for a vectorial 
operator that describes the lattice. The strength of the lattice is given by A while lo 2 is one 
of the three Bianchi VIIo invariant one-forms 

uji = dx\ , uj 2 + iuJ3 = e vpxi (dx 2 + idx 3 ) . (2) 

We see that p is the pitch of the helical structure of B^°\ The helix rotates in the x 2 -x 3 

plane as a function of the x\ direction. This helical lattice breaks translation invariance in 

the x\ direction while remaining invariant under the non-abelian Bianchi VIIo symmetry 

algebra. Our discussions of metallic and insulating behavior in the following will always 

refer to currents in the x\ direction. Currents in the x 2 and xz directions are not relaxed. 

This extreme anisotropy is a consequence of the symmetries of our microscopic lattice. We 

1 David Vegh has independently studied the conductivity in homogeneous systems that are not translation 
invariant. 
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adopt this lattice for technical convenience, but it may also be of intrinsic interest in relation 
to smectic metals [17]. An introduction to the geometry of the Bianchi VIIo algebra with 
application to chiral nematics may be found in [18]. Such helices have previously been 
found to emerge spontaneously in holographic studies including [19, 20, 21]. Here we are 
explicitly sourcing the helical structure as a simplified way to break translation invariance. 
An RG flow generated by a helical current source was previously considered in [22]. 

To extend the sources (1) into the holographic radial bulk dimension r it is sufficient to 
consider the bulk ansatz 

The bulk 

A = a(r)dt , B = w(r)uj 2 , model 

ds 2 = _t/( r )^2 + + e 2fi(r) w 2 + e 2t*(r) w 2 + e **{r)Jl (3) 

To determine the six radial functions a, w, Vi, U we need an action for the bulk theory. We 
consider the following simple Einstein-Maxwell-Proca action 

S = J d 5 x^ (r + 12- - A F ab F ab - - A W ab W ab ~ ^BaB") — ^ J B A F AW . (4) 

Here F = dA and W = dB are the field strengths. The final Chern-Simons term is not 
logically essential for our purposes, but will turn out to facilitate finding a metal-insulator 
transition. We have set the Newton constant, AdS radius and Maxwell couplings to unity. 
In general all these constants can be scaled into the fields. The wedge product is normalized 
so that the Chern-Simons term on the ansatz (3) is \pKW 2 a' . The resulting background 
equations of motion for {a,w, Vi, U} are listed in the supplementary material. 

Near the asymptotic UV boundary at r — > oo we require that the spacetime tend to 
AdS§ together with the sources (1). That is, to leading order at large r 

U = r 2 , Vi = log r , a = fi, w = A . (5) 

More important however is the far interior of the spacetime, which will describe the universal 
low energy physics [8]. Underlying the metal-insulator transition that we will find is the 
fact that there are different possible far interior geometries at zero temperature. These 

Metallic 

correspond to different solutions to the Einstein-Maxwell-Proca action. The metallic state 

GEOMETRY 

at zero temperature will be described by the well known AdSi x M 3 solution to the equations 
of motion as r — > 0, together with specific deformations: 

U = 12r 2 (l + uir 5 ), v i = v (l + v il r s ), a = 2>/6r(l + a ir s ) , w = wxr s . (6) 

The exponents are determined by the equations of motion and come in pairs 5± with 5+ + 
5- = —1. Each larger value S + gives the scaling dimension of a coupling describing a 
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deformation of this semi locally critical IR. We find two marginal operators with <5+ = 0. 
These correspond to rescaling x\ and {x2,X3}, respectively. The three remaining operators 
are found to have 

*+ = {l, ~\ + + m2 ) + T2P 2e ~ 2Vo - 7eP £ ' V ° ' -\ + \j\ + b 2e - 2Vo } ■ (7) 

If 

(2V6 k - pe- v °^J pe~ v ° <m 2 , (8) 

we can see that all three of these modes are irrelevant, with 5+ > 0. These marginal and 
irrelevant modes must be turned on and their coefficients {u±, vn,a±, w\\ tuned in order to 
match onto the near boundary behavior (5) upon integrating the equations of motion out to 
large r. In cases where (8) is satisfied, the interior solution is translation invariant, because 
w = as r — > and the deformation is turned off, and therefore the helical structure 
is irrelevant. The IR irrelevance of translation invariance breaking operators implies that 
phases with these interior geometries will be metals with a coherent Drude peak [12]. 

There is no guarantee, however, that the translation invariant solution with irrelevant 
deformations (6) is indeed the zero temperature IR geometry. The simplest way in which 
that geometry can fail to be the true IR is if the mass of the Proca field B and the Chern- 
Simons coupling are such that the inequality (8) is violated for some v Q , at fixed pitch p. 2 
This violation leads to (5+ < for one of the operators about the IR fixed point. The IR 
datum v Q can be tuned by varying UV data such as the chemical potential. Crossing the 
bound (8) will trigger a quantum phase transition as the operator becomes relevant and 
drives an RG flow away from the AdSz x M 3 fixed point. This transition could be first 
order or continuous, depending on whether the true IR fixed point is close or not to the 
RG-unstable AdS2 x M 3 fixed point when the operator first becomes relevant. This scenario 
is illustrated in figure 2 below. 

If the dimension 5+ becomes complex in (7), this indicates that the interior geometry 
has become dynamically unstable. These are the inhomogeneous instabilities investigated 
in e.g. [19, 21]. In this work we are interested in the effects of an explicit helical lattice 
rather than spontaneously generated helices. We will therefore make sure to stay at values 
of m 2 , k and p such that 5 + about the metallic AdS-z x M 3 background is never complex for 
any v Q . For the massless case that we will focus on below, the condition for stability is seen 
to be from (7) that \k\ < l/y/2. 



2 We should require that the asymptotically AdSs Breitenlohner-Freedman bound is satisfied: — 1 < m 2 . 
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Figure 2: Two renormalization group flow scenarios that arise in our theories, mediating 
quantum phase transitions between metallic and insulating phases. In the left plot, the 
phase transition is mediated by an unstable fixed point, which has relevant operators (that 
may in addition have complex scaling dimensions). In the right plot, the phase transition 
occurs when the metallic fixed point itself develops a relevant deformation. One of these 
two possibilities must occur for the metallic phase to become unstable. 

A second way that the IR can be distinct from the translation invariant solution (6) is 
if there is a bifurcation in the RG flow as a function of the UV parameters of the system 
such as the lattice strength. This is the scenario that was found to mediate holographic 
fractionalization quantum phase transitions in [23] . The bifurcation fixed point should have 
a relevant deformation and thus be unstable under RG flow. If the fixed point mediating 
the bifurcation is furthermore dynamically unstable (that is, with deformations that have 
complex scaling dimension), then the metal-insulator quantum phase transition will be first 
order, otherwise it will be continuous [23]. The required unstable fixed points are similar to 
those constructed in [22]; we are able to charge up those solutions using the Chern-Simons 
term. The Chern-Simons coupling allows the 'magnetic' B field to generate an electric A 
field with a minimal back reaction on the solution. It is simple to see that our theory admits 
the scaling solution 

U = u r 2 , vi = via , e V2 = e V2 °r a , e vs = e V3 °r a , a = a r , w = w r a . (9) 

Plugging this ansatz into the equations of motion, the exponent a and the various prefactors 
can be determined numerically by solving algebraic equations, and depend upon the pitch 
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geometry 
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p, mass squared m 2 and Chern-Simons coupling k. For the case m = that we will focus on 
in the numerics below, we find that for all p, the bifurcating solution exists for < 0.57. 

To determine the stability of this fixed point we need to perturb the functions {U, Vi, a, w} 
in (9) by a power law r s , analogously to what we did in (6) above. The operators captured 
by our bulk theory are found to now correspond to two irrelevant operators with 5+ > 0, 
one marginal operator with 5 + = and two relevant operators. This fixed point is therefore 
unstable under RG flow, as anticipated above. Depending on the regime of parameters, 
one of the relevant operators can have a complex scaling exponent. In the massless case 
m = 0, the operator dimensions are found to be independent of the pitch p, with the com- 
plex scaling exponent present for \n\ < 0.3. Thus depending on the Chern-Simons coupling, 
a metal-insulator transition mediated by this fixed point is expected [23] to be first order 
(for < \k\ < 0.3) or continuous (for 0.3 < \k\ < 0.57). Let use repeat that the solution we 
are describing at the moment should not be understood as the generic true IR ground state, 
but rather as an unstable fixed point that mediates the transition between two different 
phases. This scenario is also illustrated in figure 2. 

For the massless case m = 0, we have found the stable IR geometry that will pertain 
in the insulating phase. Unlike the metallic IR (6) and the unstable critical point (9), the 

Insulating 

interior geometry in this case is not a scaling solution on its own, but rather describes the 

GEOMETRY 

leading order terms in an expansion as r — > 0. The background to leading order as r — > 
is found to be 

U = u a r 2 + ■■■ , e V2 = e^r" 1 / 3 + • • • , e" 2 = e V2o r 2 / 3 + ■■■ , e v:i = e^r 1 / 3 + • • • , 
a = a G r 5/3 -\ , w = w D + wir 4/3 -\ . (10) 

We have quoted the subleading term in w, as this is necessary to establish consistently the 
leading order behavior of a. We see that the translation invariance breaking 'magnetic' 
background field W does not vanish in this IR background, unlike in the metallic IR (6). 
Furthermore, while the metallic IR solution (6) has electric flux, J ffi3 *F through the horizon, 
the insulating IR background (10) and the bifurcating scaling solution (9) have no such 
flux. The electric field is generated outside the horizon via the Chern-Simons term. We are 
placing the Chern-Simons term on the right hand side of the Maxwell equation and thinking 
of it as a source for charge. Our phase transition might therefore be thought of as a non- 
translationally invariant version of the fractionalization transition discussed in [24]. There, 
a Chern-Simons term together with a magnetic field enabled a phase transition between 
solutions with the charge behind and outside the horizon. While in the models we study in 
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this paper we therefore have that the metals are fractionalized and the insulators 'cohesive' 
[25], this is not essential to our setup. We expect similar metal-insulator transitions to occur 
without the Chern-Simons term, and therefore without the possibility of cohesive phases. 
In fact, we have found strong evidence for the presence of such transitions (without the 
Chern-Simons term) using numerics at finite temperature. However, we have not been able 
to identify analytically the zero temperature IR insulating behavior in these cases. Finally, 
we should note that the solution (10) will also describe the zero temperature IR of the 
helical black holes discussed in [21]. 

The stability of the insulating IR geometry (10) under RG flow is established by per- 
turbing the background by power law modes as previously. Because the background is not 
an exact solution, this is a little more subtle than previous cases. Details may be found 
in the supplementary material. The upshot of the analysis is that the background has no 
relevant modes and is therefore a stable IR fixed point. By following the irrelevant and 
marginal modes upwards along the RG flow we can construct asymptotically AdS§ domain 
wall solutions. We now have all the ingredients in place to realize in principle both of the 
metal-insulator transition scenarios outlined in figure 2. 

From the numerically constructed domain wall solutions, one can read off the asymptotic 
constants of integration from the behavior of the solutions near the AdS§ boundary as 
r — > oo. For the case of a massless B field, these take the general form 

v , /3-Ap 2 /21ogr 
a = n+^ + ---, w = A + - ^ — + ■••, (11) 



2 e/3+p 2 A 2 /6 logr 9i + Si X 2 p 2 /24 logr 
U = r 1 , Vi = log r -\ h 



Here g\ + 52 + 93 = and si = S3 = 1, S2 = —2. The logarithms indicate the presence of 
a scaling anomaly. In writing these expansions we have implicitly fixed a scale appearing 
due to this anomaly. This anomaly is a minor inconvenience, that we treat correctly via 
holographic renormalization, and will not play a significant role in our discussion. 

The free energy density is then obtained from the renormalized on shell action to be 

^ = e - nQ-Ts. (12) Phase 

STRUCTURE 

Where the charge density Q = —2v + Kp\ 2 /2. We see that the background 'magnetic' field 
B has shifted the expression for the charge density, due to the Chern-Simons term in the 
action (4). We have also allowed ourselves here to heat up the solutions. The temperature 
T and entropy density s are defined in the usual way via the surface gravity and the area 
of the event horizon. An important observable shortly will be the temperature dependence 
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of the conductivity and resistivity. In the supplementary material we describe the standard 
numerical shooting method that allows us to construct the finite temperature black hole 
solutions to our model. Including the temperature, our solutions are described by four 
dimensionful UV parameters: temperature T, chemical potential /x, lattice strength A and 
helical pitch p. Taking dimensionless ratios, this gives us a three dimensional phase space 
to explore. 

In the remainder we will focus on the phase structure at low and zero temperature. 
Firstly we present thermodynamic evidence for a continuous phase transition between zero 
temperature metallic and insulating phases described by the IR geometries (6) and (10) 
respectively. These phases are not distinguished by any symmetry breaking condensate. The 
absence of an order parameter is no obstacle in our approach, as we have a precise dynamical 
characterization of the two phases. In the metallic phase, operators breaking translation 
invariance are irrelevant at low energies, while in the insulating phase the breaking of 
translational invariance persists to the lowest energy scales. This distinction between metals 
and insulators according to whether "generalized umklapp" [12] processes are irrelevant or 
relevant at low energies has previously arisen in a soluble one dimensional model [27]. Below 
we will compute the optical conductivities and d.c. resistivities of these two phases. 

Figure 3 shows the free energy as a function of the pitch p at a fixed lattice strength and 
chemical potential. The free energy is computed at a fixed low temperature over the whole 
range and at zero temperature for some of the range. Issues with numerical accuracy close 
to the critical point have prevented us from scanning the whole range at zero temperature. 
The low temperature result is seen to very closely approximate the zero temperature result 
where the results overlap. Figure 3 shows that at small pitch the zero temperature limit of 
the system has the insulating near horizon geometry, while at large pitch the IR geometry is 
the metallic AdS2 x M 3 solution. There is therefore necessarily a quantum phase transition 
at intermediate helical pitch. We exhibit the phase transition in the figure by plotting the 
entropy density as a function of the pitch, zoomed in near the putative critical point. The 
entropy density is a useful quantity to pick out the location of the transition because while 
the AdSz x M 3 metallic geometry leads to a nonzero entropy density at T = 0, discussed 
further below, the insulating geometry (10) is easily seen to correspond to a low temperature 
entropy density that vanishes as s ~ T 2 / 3 . 

Figure 3 shows one of many cuts across the phase diagram that contain the physics of 
interest. We postpone a full exploration of the phase diagram, in particular the p — A plane 
at zero temperature, to future studies. The fact that we have chosen the parameters in the 
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Figure 3: The left hand plot shows the free energy as a function of pitch. The solid lines 
are results for T = 0, while the dashed line shows the thermodynamically preferred phase 
at T/fi = 1CP 3 . The solid line to the left of the plot corresponds to domain walls with an 
insulating IR geometry, while the solid line to the right of the plot corresponds to domain 
walls with metallic IR geometry. There is a continuous quantum phase transition between 
these two phases. The right hand plot shows the entropy density as a function of the pitch, 
for values of the pitch close to the critical value p c / fi ~ 2.1. As the temperature is lowered, 
temperatures shown from top to bottom are T / fi = {10 -3 , 10 -4 , 10~ 5 }, this curve develops 
a step function at the critical value. Both plots have m 2 = 0, K = 1/V2 and X/fj, = 3/2. 

Lagrangian to be m 2 = and k = 1/V2 ~ 0.7 > 0.57 implies that the bifurcating solution 
(9) does not exist in this instance. Therefore the quantum phase transition will realize the 
right hand case in figure 2, and is driven by an operator in the metallic phase becoming 
relevant. The phase transition appears to be infinite order in this case. It is similar but not 
identical to the quantum BKT transitions discussed in [28] . The quantum BKT transitions 
occur as fixed points annihilate at the BF bound. Here instead the dimension of an operator 
at a fixed point is tuned from irrelevance to relevance (but stays well away from the BF 
bound) . The transition is infinite order because at the critical itself the operator is marginal 
and therefore the flow away from the unstable point is maximally slow. 

To actually determine if a given solution describes a metal or an insulator we should 
of course compute the electrical conductivity. The frequency dependent conductivity is 

Optical 

obtained by perturbing the background by 

CONDUCTIV- 

5 A = e'^A^oJi , 5B = e - iwt B{r)oj 3 , 5ds 2 = e~ iut (c(r)dt®u}i+D(r)oj2®oj3) ■ (13) rrv 
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The ordinary differential equations for {^4, B, C, D} are easy to obtain but somewhat messy. 
To compute the conductivity we must solve these equations subject to the usual infalling 
boundary conditions at the horizon [26]. Near the boundary r — > oo, the Maxwell field will 
behave like 

A m + !A(°) w 2 logr , , 

A = A® + 2 — 2 — + • • • . (14) 

The logarithm is again indicative of a scaling anomaly in the boundary theory. The (holo- 
graphically renormalized) electrical conductivity is then, see e.g. [29] , 

(15, 



icoA(°) 2 

This ratio must be taken with all sources for the other modes {B, C, D} turned off. 

Our results for the real - dissipative - part of the optical conductivity are shown in 
figure 4 below. We show the frequency dependence of the conductivity for several temper- 
atures in both a metallic and an insulating phase. As we anticipated, the metallic phase is 
characterized by a beautifully conventional Drude peak whose height increases as the tem- 
perature is lowered. In the insulating phase, also as anticipated, there is no Drude peak at 
low temperatures. At intermediately low temperatures we have an incoherent metal, while 
at the lowest temperatures a soft gap opens as the d.c. conductivity does to zero. Thus 
we have realized the three phases sketched in figure 1 above. These three behaviors of the 
conductivity are widely observed in optical studies of metal-insulator transitions [3, 4]. 

The power law vanishing of the optical conductivity at zero temperature in the insulating 
phase may be computed analytically. In the supplementary material we show, using the 
matching technique of [30] , that the dissipative conductivity behaves as 

Recr ~ w 4/3 . (w->0, T = 0, insulating) (16) 

This power law agrees with our numerics. Unlike previously found power law optical conduc- 
tivities in finite density holography, this power law is not accompanied by a delta function 
at co = 0. This fact is checked numerically by plotting the imaginary part of the optical con- 
ductivity at small frequencies. We are describing an insulating phase. A power law rather 
than exponential gap in the optical conductivity is presumably related to the persistence 
of gapless degrees of freedom in the dual theory. These degrees of freedom are manifested 
in our observation above that the low temperature entropy density scales like a power law, 
s ~ T 2 / 3 , in the insulating phase. We discuss this fact further at the end of the paper. 

The second notable feature of the insulating plots in figure 4 is that the spectral weight 
from the translation-invariance delta function has been transferred both up to 'interband' 
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Figure 4: Optical conductivity in metallic (left) and insulating (right) phases. The three 
temperatures shown in each plot are T/fi = 0.011, 0.05, 0.12 (top to bottom in the metallic 
case, bottom to top in the insulating case). The metallic phase has p/jx = 2.5 and the 
insulating phase, p/p = l/y/2. Both have m 2 = 0, k = l/y/2 and X/p = 3/2. The metallic 
phase exhibits a Drude peak while in the insulating phase low energy spectral weight is 
transferred into a 'mid-infrared' peak and to interband energy scales. Note the different 
range of the axes on the two plots. 

scales of u ~ p - more specifically, the spectral weight redistribution occurs among fre- 
quencies U) < 3fj, - and also into a 'mid-infrared peak'. Such major spectral weight 
transfer, including in some cases to intermediate frequency peaks, is an experimental hall- 
mark of correlation-driven metal-insulator transitions in the cuprates and other materials 
[3, 4, 10, 11]. The intermediate peak associates an energy scale to the insulating phase. In 
cases where the metal-instulator transition is continuous, a more comprehensive study of 
conductivites across the phase diagram should see this energy scale collapse at the critical 
point. The metallic phase is featureless at intermediate energy scales. 

Figure 5 shows the d.c. resistivity p = 1/a as a function of temperature for two insulating 
and three metallic systems. At low temperatures the resistivity increases as a function of 

D.C. 

temperature in the metallic phase and decreases in the insulating phase. Indeed these are 

RESISTIVITY 

the defining features of metals and insulators, respectively! 

The low temperature behavior of the plots in figure 5 can be understood analytically. In 
the insulating phase we expect the resistivity to diverge as p ~ j 1 - 4 / 3 ) corresponding to the 
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Figure 5: Log-log plot of the d.c. resistivity as a function of temperature. From top to 
bottom, the upper (blue) lines are insulating phases with p/p, = 1,1 /y/2 and the lower 
(red) lines are metallic phases with p/p = 2.2,2.4,2.5. As before, m 2 = 0, k = l/v2 and 
X/p = 3/2. The resistivity po is an arbitrary scale. 

scaling found in (16). In the metallic phase, the analysis of [12] implies that the resistivity 
will go to zero like p ~ T 2 ^ p \ where A(p) is the smaller of the irrelevant exponents 5+ in 
(7). This is the least irrelevant IR operator responsible for breaking translation invariance. 

The UV completion of our boundary theory is a 3+1 dimensional CFT. Therefore, at 
sufficiently high temperatures dimensional analysis requires that the resistivity decrease 
like 1/T in all phases. This feature can be seen in our numerical results, not shown here, 
but is quite independent of the universal IR dynamics that is our primary interest. It 
has the significant consequence, however, that in the metallic phases there is a maximal 
resistivity at intermediate temperatures. While this point remains to be systematically 
investigated, there does not seem to be a bound on how large the resistivity can get while 
the system remains metallic (that is, while the resistivity is still growing with temperature). 
We therefore obtain, to our knowledge for the first time, a controlled demonstration of the 
theoretical possibility of bad metals once a quasiparticle picture of transport is abandoned 
[7, 4]. The tight connection we find between bad metals and metal-insulator transitions 
resonates strongly with the experimental facts [4, 5, 6]. 

We have demonstrated that interaction-driven metal-insulator transitions can be realized 
with relative ease in a holographic framework. This fact opens many exciting possibilities 
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for connecting holography with some of the deeper challenges in contemporary condensed 
matter physics. Within our model, it remains to complete our exploration of the phase 
diagram and exhibit the other class of bifurcating quantum phase transitions, perhaps as a 
function of the lattice strength A. Looking beyond our model, let us mention two natural, 
concrete extensions of this work. 

Our metallic phase had a zero temperature IR geometry described by AdS2 x M 3 . It 
is well known that such phases exhibit a ground state entropy density and much of their 

Extensions 

interesting 'semi locally critical' [31] phenomenology may therefore be an artifact of the 

OF OUR 

large N limit of holography. It has recently been emphasized, however, that semi local 

MODEL 

criticality can also easily arise in holographic models without a ground state entropy density 
[32, 33, 34]. In particular, these models will continue to exhibit power law resistivities as a 
function of temperature due to lattice scattering in the metallic phase [12, 34]. The essential 
physics we have discussed can be adapted to these models and it is of interest to see how 
this affects the onset of insulating behavior. 

An interesting if exotic feature of our model is that the insulating phases exhibited a 
soft gap in the optical conductivity, with a power law suppressed spectral weight (16). A 
stronger suppression of spectral weight could be achieved if the full theory were gapped in the 
insulating phase, for instance by confining the SU(N) gauge fields of the dual field theory. 
The suggestion that holographic confining phases should describe Mott insulators was made 
in e.g. [35, 36, 37]. These papers did not break translation invariance and therefore did not 
exhibit insulating behavior. However, it would be natural to adapt our analysis to a system 
in which the onset of insulating behavior was tied to a confining phase transition. Indeed 
it seems plausible that gapped theories with broken translational invariance will always 
be insulating. Related to this point, it is also of interest to understand to what extent 
insulating behavior can be decoupled from the appearance of cohesive charge carriers [25]. 

The insulators obtained in this paper do not appear to be Mott-like in the sense that 
the insulating phases are not tied to charge commensurability. This raises the possibility 
that we have described a qualitatively new type of strongly correlated gapless insulator. 
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Supplementary material 

The equations of motion for the background 

The equations of motion following from the bulk action (4) on the metric and vector field 
ansatz (3) are 

= a" + a' (v[ + v 2 + v' 3 )+ P k e^ vl+V2+v ^ ww' , (17) 
= w" + w' (^ + v[-v' 2 + v' 3 ^J - (m 2 + pV 2 ^-^) - p n a! e'^^) ^ , 
= a' 2 - e~ 2v2 Uw' 2 + 2U\v[+v' 2 + v' 3 ) + 4U{v[v' 2 + v' l v' 3 + v' 2 v' 3 ) 
-24 + Ap 2 e- 2vi smh 2 (v 2 - v 3 ) + w 2 (m 2 e~ 2v2 + p 2 e " 2 K+^ , 

o = 4K + <) + ^ + e - 2 -^ 2 -^ + 4(^K + 4) + ,; 2 + ^ + ^) 

-— + (m 2 e~ 2v2 - P 2 e~ 2{ - V1+V3 ^j — - (2e~ 2vi + e - 2 («i-»2+i>3) _ 3 e -2(vi+u2-u 3 )j ?! ; 

= ^v' 3 ^v';) + 2 ^-e- 2 ^w' 2 -^+4(^(v[ + v' 3 ) + v[ 2 + vi + vW^ 

_ ^ m 2 e" 21 ' 2 -^e" 2 ^ 1 ^ 3 ^ — - ^e -2 " 1 - 3 e - 2 (w-^+«3) + e -2(yi+v2-v 3 )^ ^ ; 
= 4K + 40 + ^ + e- 2 ^ 2 -^ + 4^(^ + ^)+^ + ^ + ^) 

+ ^ m 2 e " 2t, 2 _ p 2 e -2( t , 1 +^)^ ^ _ ^_ 2e - 2i 'i + e -2(vi-vz+v 3 ) + e -2{ Vl +v 2 -v 3 )^ P_ _ 



There is a further second order equation of motion that follows from the above equations. 
Note in particular that the third equation above is a first order equation. 

Finite temperature domain walls 

We demand the existence of a regular black hole horizon in the geometry at the radius 
r = r + at which U(r+) = and the temperature T = U'(r+)/ (47r) is finite. We are led to 
the analytic expansion of our functions 

U =l7 ( i) (r - r+) + ■ ■ ■ , 

w =w (0 ) + W(i) (r - r+) H , 

a =a (0 ) (r - r+) + ■ ■ ■ , 

Vi =Vi( ) +«i(i) (r-r + ) H , i = 1,2,3, (18) 

where dots denote higher order terms in r — r+. All of the expansion coefficients in (18) 
are determined in terms of the six constants of integration {tf(o) ; O(o)> ^(o)) w i(i)} v ^ a the 
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equations of motion. In general we will need 11 constants of integration for a unique solution 
of the differential system (17). On the other hand, in our asymptotic expansion (11) we 
find another five independent constants of integration, for given values of the sources at 
the boundary, that can be taken to be {u, (3, e, 51,52}- This gives a total of 11 constants of 
integration for a given r + , which we can keep unfixed in order to move through solutions of 
different temperature. 

After fixing r + and the deformation parameters {//, X,p} we now have a two point 
boundary value problem to solve. This will give at most a discrete set of solutions for the 
11 constants {^(o), O(o)> v i(o)> v i Pi e > Pi> #2}- We have solved for these constants using 
a shooting method. 

Domain wall solutions for insulating and metallic phases 

The zero temperature domain walls are constructed in the same way as the finite tempera- 
ture solutions just discussed. The only difference is the form of the near-horizon expansions. 



The insulating IR 

In expanding our solutions in the IR we find it is necessary to use a slightly different 
coordinate choice to that presented in the main text. Relative to (10) we let r — > r + r+. 
The expansion we use for our fields close to the insulating IR solution then reads 

U = L 2 (r - r+) 2 j 1 + (r - r+) i/3 + £ t£ (r - r+) Si + ■ ■ 

\ 3 j 

Vi =v m + Pi ln(r - r+ ) + Y, df {r - r+) j/3 + ^ cj. (r - r+) 5 * + 

3 3 

a =(r - r+) d% (r - r+)^ 3 + ^ cj (r - r+) Si + • ■ 

\ 3 3 

w =w (0) + (r - r + ) 2 / 3 dj (r - r + ) j/3 + £ cj, (r - r + ) 5i + • • • J , (19) 

where L = Sy^, pi = -1/3, p 2 = 2/3, p 3 = 1/3, e 1 ^) = 2e v ^) +v ^/p and w (0) = 
A^e 2v ^) +V ^ /p 2 . In the above expansion the constants d represent the background ge- 
ometry that governs the near horizon limit and are uniquely fixed by solving the equations 
of motion 3 . The c's are constants of integration, of which the three constants {c^c 2 , , c 3 } 



3 The expansion (19) will be true in general. However, there is a discreet set of values for the CS coupling 
k for which our expansion (19) for the background will be modified by logarithmic terms. 
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can be chosen at will. These are the deformations parameters corresponding to the three 
irrelevant operators of the IR theory with 

<5i = I (-5 + \/l + 120 K 2 ) , 5 2 = * (-5 + v 7 !^) , <5 3 = * (-5 + v 7 !^) . (20) 

Finally, the constants {ui(o) 5 w 2(o)> r +} that we have included at leading order in the expan- 
sion correspond to marginal modes. We therefore have six constants of integration with, as 
previously, another five constants from our UV expansion (after fixing all deformations and 
the pitch of the helix) . We thus expect to have a unique solution. 

The metallic IR 

In the metallic phase we again employ a shift of coordinates r — > r + r + relative to those 
used in the main text in (6). In the case where we have AdS% x M 3 as the IR limit of our 
zero temperature solutions the near horizon expansion reads 

U =12 (r - r+) 2 + Yl <% ( r " r +) &l+2 + ■ • • > 

i 

Vi =V m + Yl ^ ( r ~ r +^' + " ' ' 
i 

a =2VQ (r -r + ) + Y 4 {r ~ r+f l+1 , 

i 

^=E C ^ (r~r + f , (21) 

i 

where the constants {c^c^Cg 2 } can be chosen freely and correspond to deformation pa- 
rameters of the IR theory with 

5 1 = 1, < j 2 = _i+ v /i + ^p2 c -2«i_^p C -« lj 5 3 = 4 + ^ + |e- 2 V. (22) 

These are the exponents quoted in (7) in the main text, here with the mass m? = 0. The 
constants {vu \, v 2 m) = f3(o)> r +} which appear in the expansion represent marginal modes 
of the AdS2 x M 3 solution. Once again there are six free constants in the IR, leading to a 
unique solution for given deformation parameters in the UV. 

Holographic renormalization of action 

As usual, our bulk action (4) needs to be supplemented by appropriate boundary countert- 
erms. In particular we need to add 



Sbdy = d x v 73 ! 



2 K - 6 + 1 In r {F mn F mn + G mn G mn ) + ■■■ 



(23) 
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Here K = 7 mn V m n n is the trace of the extrinsic curvature of the boundary, n m is an 
outward pointing unit normal vector and 7 is the determinant of the induced metric at 
infinity. The logarithmic terms cancel the divergence related to the trace anomaly in four 
dimensions. The ellipsis indicate a Ricci scalar term of the four dimensional boundary 
metric, which will not contribute given our asymptotics (11). 

Analytic results for the optical conductivity 

At the lowest frequencies, the dissipative part of the conductivity (15) can be computed 
purely from near-horizon IR data. The general matching formula derived in [30] says that 
as co — > 

R^( w ) = £d'^M. (24) 

Here the sum is over all operators 0/ in the IR theory described by the near horizon 
geometry. The operators are taken to diagonalize the IR scaling operator. That is, they 
are decoupled at a quadratic level in the bulk. The d 1 are real coefficients that are nonzero 
if the corresponding operator has an overlap with the current operator J in the UV. The 
G^(ui) are the retarded Green's functions of these operators in the IR theory. We must 
therefore proceed to obtain the Green's functions for the IR insulating geometry (10). 

The equations of motion for the frequency dependent perturbation (13) about the insu- 
lating IR background (10) are fairly simple. Keeping the lowest powers in r — > one finds 
immediately a decoupled equation for the Maxwell perturbation A(r) 

A further mode is decoupled by defining 



E = *i£r» (26) 

which is then found to satisfy 

w/ 5(288 r 2 -5 u 2 ) „ 

E" - -i — '-E = . 27) 

324 r 4 v 1 

These two equations are easily solved in terms of Hankel functions. Imposing infalling 
boundary conditions at the horizon, 



(28) 
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Expanding for r — > oo and taking the ratio of the normalizable by the non-normalizable 
mode, as required in the matching procedure [30], one finds 

Im g%(u) ~ lj 7 / 3 , Im Q§(uj) ~ w 13/3 . (29) 

There is a third physical mode in the IR that is more difficult to decouple. However, because 
the IR equations of motion only depend on the ratio co/r, it is possible to extract the power 
of uj by solving the equations of motion with u = 0. This is equivalent to finding the IR 
dimension of the operator. With oj = the remaining IR equation is 

2 40 

B" + -B' - -^B = . (30) 

This has solution B = cir 5 / 3 + C2r~ 8 / 3 . Dimensional analysis then implies that this operator 
will also lead to 

lmg§(oj) ~w 13 / 3 . (31) 

From (29) and (31) we see that the most singular of the three IR operators is A. Thus we 
conclude from (24) that the optical conductivity at small frequencies will scale as 

Rea(w)~w 4/3 . (32) 

This is the result quoted in the main text in (16) and observed in our numerical results. 
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